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Abstract 

Some simple guidelines based on the accuracy in deter- 
mining a satellite formation's semi-major axis differences 
are useful in making preliminary assessments of the nav- 
igation accuracy needed to support such missions. These 
guidelines are valid for any elliptical orbit, regardless of 
eccentricity. Although maneuvers required for formation 
establishment, reconfiguration, and station-keeping require 
accurate prediction of the state estimate to the maneuver 
time, and hence are directly affected by errors in all the or- 
bital elements, experience has shown that determination of 
orbit plane orientation and orbit shape to acceptable levels 
is less challenging than the determination of orbital period 
or semi-major axis. Furthermore, any differences among 
the member’s semi-major axes are undesirable for a satel- 
lite formation, since it will lead to differential along-track 
drift due to period differences. Since inevitable navigation 
errors prevent these differences from ever being zero, one 
may use the guidelines this paper presents to determine 
how much drift will result from a given relative naviga- 
tion accuracy, or conversely what navigation accuracy is 
required to limit drift to a given rate. Since the guidelines 
do not account for non-two-body perturbations, they may 
be viewed as useful preliminary design tools, rather than 
as the basis for mission navigation requirements, which 
should be based on detailed analysis of the mission con- 
figuration, including all relevant sources of uncertainty. 

Introduction 

One of the most significant differences between many for- 
mation flying satellite missions that are currently of inter- 
est, and the intentionally close approaches that past mis- 
sions have performed - e.g. rendezvous, docking, and prox- 
imity operations of the Space Shuttle - is the need for 
long-term and efficient maintenance of relatively close for- 
mations. As of 2002, separations of a few hundred kilome- 
ters have been achieved for long-duration formations, and 
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in the next few years separations in the range of 100 m to 
10 km are expected to be attempted. Ref. 1 describes some 
relevant aspects of one such recent mission. 

Many perturbations affect the accuracy of maneuvers re- 
quired to maintain these formations, but in principle the 
maneuver planning process can accommodate all known 
perturbations. However, the navigation errors at the time of 
the maneuver computation will form a lower bound on the 
accuracy of the maneuvers, since even a perfectly executed 
maneuver will ‘flock in” the navigation errors. Indeed, 
Ref. 2 suggests that velocity uncertainty may be the lim- 
iting technology for formation flying. 

Although the errors in all the states will affect the maneu- 
ver, as Ref. 3 discusses, semi-major axis error is the most 
difficult to estimate. An essential point of Ref. 3 is that 
(for circular orbits) semi-major axis uncertainty depends on 
three quantities: the radial position error, the along-track 
velocity error, and the correlation between these errors, 
which arises due to the the conservation of energy. The 
most notable consequence of semi-major axis differences 
for a formation of satellites is relative drift in the direc- 
tion of the orbital path, which Figure 2 illustrates. For 
the relatively high-thrust, short burning propulsion systems 
flown on current missions, the result of this drift will be 
more frequent stationkeeping maneuvers. For “continuous” 
low-thrust systems proposed for many upcoming missions, 
this drift will be random error that must be counteracted, 
limited by the accuracy of the navigation feedback signal. 
Figure 1 illustrates the effect of the various contributors 
to along-track drift in a circular orbit, based on the rela- 
tionships of Ref. 3. To interpret Figure 1, consider a nav- 
igation system that can produce radial accuracy of 10 cm, 
and speed accuracy of 0.1 mm/sec. Figure 1 indicates that 
the corresponding along-track drift may be on the order of 
1 m/ orbit, and could be as poor as about 3 m/orbit, if the 
correlation between radius and speed, p rvy is poor. 

Ref. 4 describes a necessary condition for obtaining a 
no-drift solution in eccentric orbits in a deterministic set- 
ting. Ref. 1 describes requirements on relative semi-major 
axis knowledge and control for relative drift of the GRACE 
mission due to differential drag. The contribution of the 
development below is to describe and quantify the relative 
drift issue interns of the (non-deterministic) navigation er- 
tots for elliptical orbits, generalizing the results of Ref. 3. 
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Fig. 1 Along-track drift due to semi-major axis error for a typical low earth orbit 3 Each family of contours is based on a 
constant semi-major axis error, resulting from various combinations of radial position error, along-track velocity error, and the 
correlation between these. 


Relative motion due to 8a=1 1 m for e=0.8 over 15 orbits 



Tangential Separation [km] 


Fig. 2 Along-track drift due to semi-major axis differences 
between two spacecraft in a highly elliptical orbit 

The paper concludes with examples in low circular and 
high apogee elliptical Earth orbit missions. 

Development 

For every complete revolution in any two-body elliptical or- 
bit, an error in knowledge of the orbit period, 5T P> results in 
an along-track error growth proportional to the speed (ve- 
locity magnitude), v, 

( 1 ) 


where /$ is the initial true anomaly. It is not hard to show 
(viz. Ref. 3) that 



where a is the semi-major axis and fi is the gravitational 
constant, GM, so that 

Ss(fi 4 27 r) = -3 7tv(/< 4- 27t) (3) 

At periapse, 

w(0)=i;(0,2ir,...) = ^yiH (4) 

where e is the eccentricity, and so the in-track growth per 
orbit at periapse is 

$$(2^) = . .) = — 37r (5) 

V 1 — e 

while at apoapse, 



and hence the in-track growth per orbit at apoapse is 

&$(37r) = ds(37r,57r 5 . . .) = —%n\\ — -5a (7) 


6s(fi 4- 2 tt) = -v(/i 4 2t t)5T p 


2 OF 7 




Since according to Eq. (5) or (7) the along-track error 
growth per orbit is linear in the semi-major axis error, their 
error distributions share this linear relationship. For exam- 
ple their standard deviations are related by 

<7 S (7T) = 3^Y / Y“j^ Cr a (8) 

Ref. 3 gives the following relation, valid for circular orbits, 
among semi-major axis error, radius error, speed error, and 
the correlation between radius and speed errors,. p TV : 


(T 0 = 2 J <T? + —prv&rCTv + (9) 

V n n z 


any point in the orbit, the following relationships are useful: 
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where n is the mean motion. Note that the velocity terms 
in Eq. (9) are with respect to an inertial frame. Relative to 
a frame rotating at the (constant) orbital angular velocity, 
such as Hill’s frame, 


/ 4 i 

<T a = 2W 4<J 2 -f -Pry<7r&y 4 —$0% (10) 

where y in Eq. (10) refers to the along-track component of 
velocity relative to the rotating frame. 

The general form of the semi-major axis variance, also 
given in Ref. 3, is 

= A(x ref )P x A(x ref ) T ( 11 ) 

where the state is x T = [r T , v T ], r is the position, v is the 
velocity, P x is the state error covariance, x re f is given, and 


where p = \/l — e 2 . From the definition of Eq. (12), 5a = 
A(x re /)dx, which may be written in terms of radial and 
along-track components, x and y, respectively, as 



Using Eqs. (15) - (20) in Eq. (21), it is not hard to show' 
that 


5a 


= 2/- 
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Taking the expectation of the square of Eq. (22), and as- 
suming a zero mean, results in the following expression for 
the semi-major axis variance at any point in an elliptical 
orbit: 


A(x) = 2a 2 



( 12 ) 


Applying Eq. (11) at apoapse, and using Eq. (7), it is not 
hard to show that Eq. (9) may be generalized to 
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or similarly at periapse to 
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To get an expression for the semi-major axis variance at 
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If the correlations between radius and radial velocity, 
Prv x , and between radial velocity and along-track velocity, 
Pv x v y > can be assumed to be insignificant, then the follow- 
ing simpler expression results: 


<r 2 _ (1 + ecos/) 4 2 
4 “ 77 s <Tr 

+ jAy [( esin /) 2 < + (1 +ecos/) 2 ff 2 v j 

2(1 4 ecos/) 3 
' C Prv y &T &Vy 


(24) 
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e * 0.8, a = 42096 km, a = 10 m, c « 10 mm/s, o * 1 mm/s 

r vy vx 



Fig. 3 Semi-major axis uncertainty as a function of true 
anomaly for a particular high Earth orbit example. 


Figure 3 illustrates Eq. (24) for a particular example. Note 
that different relationships among a r , a Vy , and a Vx than 
Figure 3 shows produce quite different results. 

In most satellite formations, it is desirable for the satel- 
lites to minimize their semi-major axis differences so as to 
minimize relative along-track drift. The general form of 
the variance of the relative semi-major axis between any 
two satellties i and j> also given in Ref. 3, is 

= A(Xi)PiA(Xi) T + A(Xj)PjA(Xj) T 

- A(x.i)P ij A(x j ) T - A(x j )P?A(x i ) T (25) 

where — PxiXj * To approximate and simplify Eq. (25), 
assume that die satellites in the formation have approxi- 
mately the same a and e, that their navigation errors have 
the same distributions, and that Pij = pijPi. Then, the 
version of Eqs. (13) - (14) and (23) - (24) corresponding 
to relative semi-major axis errors is 

tf’Aa = ^ Pij <7 a (26) 


Therefore, as in Eq. (8), the standard deviation of the in- 
track drift per orbit between any two satellites in a forma- 
tion, evaluated at each apoapse, is related to their relative 
semi-major axis standard deviation by 

<7a s (3tt) = -3^y (27) 

In a good orbital navigation filter, —1 < p rv < -.9, 
due to die dynamical constraint imposed on the estimates 
by conservation of energy*. A consequence of the high 
correlation between radius and speed errors is that speed 
error may be viewed as a dependent variable of radius error. 
Lear 5 gives the following approximation of the speed error 
required to “balance” a corresponding radius error: 


y/T^ l + je 2 + ie 4 + 


(28) 


256 


*Ref. 3 shows that one should not always assume that GPS-based or- 
bital navigation systems properly account for this constraint 


Eq. (28) gives the speed error that produces the same size 
error in semi-major axis as the given radius error. Table 1 
gives some numerical examples of Eq. (28). If Eq. (28) 


Table I Speed uncertainty, (j v (m/s], “balancing” radius un- 
certainty for a range of eccentricities, according to Eq. (28). 


e 

0.01 

0.18 

<7 r [m] 

3.1 

56 

1000 

0 

1.2e-05 

0.00022 

0.0039 

0.070 

1.2 

0.35 

2.4e-05 

0.00044 

0.0077 

0.14 

2.5 

0.7 

9.4e-05 

0.0017 

0.030 

0.53 

9.4 

0.9 

0.00074 

0.013 

0.23 

4.2 

74 


holds, then in-track drift per orbit may be viewed as a func- 
tion of the radius error only. 

If there are no relative measurements, but instead each 
satellite’s absolute state is estimated separately, e.g. from 
GPS pseudoranges, pij may be quite small. If the filter 
processes relative measurements such as cross-link ranges 
and/or GPS meausurement differences, p%j « 0.9 may be 
a reasonable assumption. Based on the assumptions that 
Pij ~ 0.9, that the approximate velocity constraint of 
Eq. (28) is valid, and that the radius and speed are well- 
correlated {prv = —.9), Figure 4(a) gives some examples 
of relative drift due to semi-major axis error for various ec- 
centricities, for Earth-orbiting formations. The left subplot 
illustrates drift rates at apogee, and the right subplot drift 
rates at perigee. From the figure, it is clear that as eccen- 
tricity increases, the formation will drift apart more slowly 
at apogee, and more quickly at perigee, for a given radius 
or semi-major axis error. 

Figure 4(b) shows similar results, instead assuming no 
correlation between the two satellites’ state estimates. Fig- 
ure 4(c) shows some additional results, in which a “fudge 
factor” of 50% is applied to Eq. (28), and zero correlations 
assumed. This is intended to capture the type of “poor” ve- 
locity estimation that Ref. 3 describes may occur with some 
GPS-based orbit determination systems. 

Finally, it may be of interest to determine the time to 
drift a given distance due to a semi-major axis error. Since 
time-to-drift depends inversely on the drift rate, a nonlinear 
mapping of the statistical moments of the drift rate per orbit 
is required to find the statistics of the time to drift a given 
distance. Let D be the deadband size, i.e. the size of the 
“control box” in which the satellite must remain, and r be 
the time to reach the deadband. Assume that the relative 
navigation errors are zero-mean and Gaussian, then due to 
their linear relationship, the relative drift rate also has a 
Gaussian probability density, 

Since 

r = g(As) = ^ (30) 

is a nonlinear function, its probability density may be found 
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In-track relative drift uncertainty per orbit, o # [mj In-track relative drift uncertainty per orbit, 0^(m] In-track relative drift uncertainty per orbit, 


Apogee 


Perigee 


in terms of Eq. 29 via 6 



Apogee Perigee 



Apogee Perigee 



c) “Poor” velocity accuracy and zero correlations. 

Fig. 4 Relative drift per orbit due to semi-major axis error 
for various eccentricities 


fy(y) = 

1 


fx(xi) 

b'&)l 


*i=p' 1 ( y) 


(31) 


where the Xi are the values of x where y — g(x). Using 
Eqs. (29) and (30) in Eq. (31) results in 



where the fact that the deadband is reached with a positive 
or negative drift rate has been used. 

Figures 5(a) and 5(b) show the probability density and 
the probability of not reaching the deadband in time r for 
various values of D / a 

Figure 5(b) can be used for finding the minimum accept- 
able value of the deadband. For example, if one wants 
there to be a 75% probability that the deadband will not 
be reached in four orbits, then D/<j& s > 5, keeping in 
mind that this analysis does not include perturbations such 
as differential drag or thrust errors. 


Examples 

Suppose there is some interest in flying a high-eccentricity, 
e -= 0.8, formation flying mission, in which the satellites 
are supposed to be 1 0 km apart at apogee. If the separa- 
tion varies by more than 10% at apogee, or if there is some 
danger of a collision, the satellites must perform short- 
duration, relatively high-thrust stationkeeping maneuvers. 
In order to maximize the science return and mission du- 
ration, such maneuvers should not occur more often than 
every four weeks. The orbit period is approximately one 
day, so the relative drift at apogee must be less than one 
kilometer per 28 orbits, or about 36 meters per orbit. Con- 
sulting Eq. (27), the relative semi-major axis error should 
therefore have a “one sigma” value of 36/7 r — 11 me- 
ters. If relative navigation system is highly correlated, i.e. 
Prv = — -9 and pij = 0, Figure 4(a) indicates that a “one 
sigma” radius error of about 55 meters at apogee should be 
sufficient to meet the 10% requirement in a “one sigma” 
sense. If the navigation system is not well-correlated, i.e. 
p rv — p^ — 0, and the relative velocity errors dominate as 
with many of the GPS systems described in Ref. 3, then 
Figure 4(c) indicates that a “one sigma” radius error of 
about 2 meters would be necessary to meet the 10% re- 
quirement. 

Note however that at perigee, the 1 1 meter relative semi- 
major axis error would produce a drift of about 310 meters 
per orbit, which may be found using Eq. (5). From Fig- 
ure 4(a), this drift corresponds to a radius error of 0.6 m 
for a well-correlated system. For a poorly correlated sys- 
tem like Figure 4(c) illustrates, the corresponding radius 
error that would be required at perigee is about 0.15 m. 
For most elliptical orbit formations, the separation will in- 
crease at perigee, so a few kilometer change in relative 
position may amount to much less than 10%. However, 
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Number of orbits until drift to control box boundary 

a) Probability density function for time to reach deadband, for var- 
ious ratios of deadband size to drift rate. 


Probability of not reaching control box boundary 



Number of orbits 

b) Probability of not reaching the deadband in given number of or- 
bits, for various ratios of deadband size to drift rate. 

Fig. 5 Drift probabilities. 

there do exist elliptical orbit formations for which separa- 
tion at perigee decreases relative to apogee, in which case 
a collision avoidance maneuver might be required much 
more often than desired unless the navigation error could 
be further reduced. 

Next, suppose a low-Earth orbiting (550 km altitude), 
leader-follower mission is proposed in which the separa- 
tion distance is 60 m, and this must be maintained within 
±20 m. The mission will use a highly accurate differen- 
tial Global Positioning System relative navigator, which 
can determine relative position to within 6 cm, and relative 
velocity to within 2 mm/sec, “one sigma,” per axis. From 
Figure 1, one can see that the corresponding drift uncer- 
tainty is about 3 m, that the velocity noise is the limiting 
error, and that the position/velocity noise combination is 



Fig. 6 Along-track control error time history- for LEO leader- 
follower formation with semi-major axis errors. 

such that the value of the correlation coefficient has a re- 
duced effect, since the dominating velocity noise places the 
position/velocity noise combination above and to the left of 
the main diagonal of the figure, which is where the curves 
split based on correlation coefficient (this is consistent with 
the findings of Ref 2). It therefore seems appropriate to 
use the poorly correlated example of Figure 4(c) to assess 
the drift rate per orbit, from which one can determine the 
drift to be about 5 m per orbit For a 20 m deadband, one 
should therefore interpolate between the D/a& s = 3 and 
D/a A s = 5 curves of Figure 5(b), from which it is possible 
to determine that there is about an 80% probability that the 
deadband will not be reached in less than about 5 orbits, 
or equivalently that there is a 20% probability that every 
5 orbits, the deadband will be reached. 

Finally, Figure 6 demonstrates the effect that relative 
navigation error has on a formation. The figure shows the 
along-track relative motion time history for one week (105 
orbits) for a leader-follower formation with a desired sepa- 
ration of 1 00 m with the 2 mm's relative velocity and 6 cm 
relative position relative navigation errors described above. 
The relative state control strategy is a minimum fuel in- 
track impulsive strategy in which each maneuver consists 
of three impulses separated by 0.5 orbits. A maneuver is 
initiated when 75% (15 in) of the deadband of 20 m is 
reached. 7 To initiate relative motion there is differential 
drag, but after the first maneuver the differential drag is set 
to zero and the only error is the relative navigation error. 
As Figure 6 shows, when there is a high drift rate due to 
the semi-major axis error the deadband is exceeded. 

Conclusion 

This paper has presented some guidelines on relative nav- 
igation that may be useful for conceptual analysis of for- 
mation flying missions. These guidelines generalize pre- 
viously reported formulae to the case of elliptical orbits. 
Since the guidelines do not account for non-two-body per- 
turbations, they may be viewed as useful preliminary de- 
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sign tools, rather than as the basis for mission navigation 
requirements, which should be based on detailed analysis 
of the mission configuration, including all relevant sources 
of uncertainty. 
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